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SUMMARY 

The  torsional  rigidity  of  solid  cylinders  of  double-wedge  section 
is  considered  tteoretioally . Minimum  energy  methods  are  used  to 
determine  close  upper  and  lower  limits  to  the  rigidity.  The  results  arc 
presented  in  graphical  form. 
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1 Introduction 


In  this  report  the  torsional  rigidity  of  solid  cylinders  of  double 
wedge  section  is  considered  theoretically.  A lower  limit  for  the 
rigidity  has  been  obtained  in  a manner  similar  to  that  used  by  Dunoan1; 
a parabolic  variation  of  the  stress  function  across  the  thickness  is 
assumed  and  the  Ritz2  method  is  then  used  in  conjunction  with  a varia- 
tional technique  to  determine  the  rigidity.  An  upper  limit  has  been 
obtained  from  the  static  analogue  of  Kelvin's  theorem^;  a linear  variation 
of  the  warping  function  across  the  thickness  is  assumed  and  a variational 
technique  then  used  to  determine  the  rigidity. 
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List  of  Symbols  (See  Figure  1) 

C = torsional  rigidity 

G = shear  modulus 

t = maximum  thickness  of  section 

c = chord  of  section 

X = fraction  of  chord  at  which  maximum  thickness  occurs 

m = Vc  ratio 
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(m1  + m2)  (m1  - mg  + p1  (l  - m1  m2)j 
(Pi  + P2)0  “ " ®|) 


(m1  + m2)  {m2  - m1  + p1  ( 1 - m1  m2)| 
(P1  + P2)(1  " 


3 Lower  and  upper  limits  for  the  torsional  rigidity 

A lower  limit  for  the  rigidity  has  been  found  in  Appendix  I on  the 
assumption  that  the  stress  function  varies  parabolically  across  the 
thickness;  the  rigidity  is  then  determined  by  the  Ritz  method  and  a 
variational  technique.  An  upper  limit  for  the  rigidity  has  been  found 
in  Appendix  II  on  the  assumption  that  the  warping  function  varies  linearly 
across  the  thickness;  the  rigidity  is  then  determined  from  the  static 
analogue  of  Kelvin's  theorem  and  a variational  technique.  It  follows 
that  the  torsional  rigidity  satisfies  the  inequality 
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These  limits  have  been  plotted  in  Figure  2 for  various  values  of  X 
up  to  Vc  = 0.3.  It  will  be  seen  that  over  the  range  considered  the 
limits  are  dose;  the  maximum  error  that  can  arise  by  taking  the  mean 
of  the  two  limits  is  less  than  1.6)%. 


Equations  (2)  and  (3)  may  be  simplified  far  the  special  cases  in 
which  the  section  becomes  a diamond  or  a triangle. 

3.1  Special  gage:  diamond  section  (X  = 0.5) 

For  a diamond  section  equations  (2)  and  (3)  reduce  to 


Jlower 
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and 

G-ct3  1 - 5 m2  - if  nA  + if  m3  'i  3 + m2 
Gupper  ' 12  L (1-m2)2 

3 .2  Special  case:  triangular  section  ( X = 0 or  1) 

For  a triangular  section  equations  (2)  and  (3)  reduce  to 
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4 Discussion  of  Results 


It  will  be  seen  from  Figure  2 that  when  the  maximum  thickness  is 
near  the  mid-chord  (i.e.  X = 0.5)  the  torsional  rigidity  is  practically 
independent  of  X , which  is  to  be  expected  from  considerations  of 
symnetiy.  For  a cylinder  for  which  yc  < 0.05  and  0.2  < X < 0.8  the 
torsional  rigidity  is  approximately  Set 3/ 12  which,  for  materials  in 
which  v = i,  is  1.5  times  the  flexural  rigidity. 

For  a given  Vc  ratio  the  lower  and  upper  limits  are  closest  when 
the  section  is  a diamond  and  are  furthest  apart  when  the  section  is  a 
triangle.  If  Vo  = 1 and  X = 0.5  (corresponding  to  the  limiting  case  of 
a square)  the  lower  and  upper  limits  are  each  in  error  by  3.6#,  and  if 
t/c  = 1/fb  and  X = 0 (corresponding  to  the  limiting  case  of  an  equila- 
lateral  triangle)  the  lower  limit  is  correct  and  the  upper  limit  in 
error  by  12.6%. 

5 Conclusions 


The  torsional  rigidity  of  solid  cylinders  of  double-wedge  section 
has  been  considered  theoretically.  Minimum  energy  methods  have  been 
used  to  det ermine  close  upper  and  lower  limits  to  the  rigidity.  The 
variation  of  the  torsional  rigidity  with  the  t/c  ratio  and  with  the 
position  at  which  the  maximum  thickness  occurs  has  been  investigated 
and  the  results  presented  in  graphical  form. 
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Additional  Symbols  used  only  In  Appendices  (See  Figure  1) 
Ox,  Oy  j 

, 01y1  'v  Cartesian  co-ordinates 

°2X2»  °&2  J 

xsy  = co-ordinates  of  centre  of  twist 

0 = torsion  stress  function 

w = warping  stress  function 

c1  = 

c2  = ( 1 -X) c 

K,  H = surfaoe  integrals 

f , g^=  functions  of  x1 

f2>  g2=  functions  of  x2 
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jffiHPiDIX  I 

Calculation  of  lower  Unit 

In  the  Rite  method  a form  for  the  stress  function  <f>  Is  chosen 
that  vanishes  on  the  boundary  of  the  section  and  which  nay  contain  a 
number  of  arbitrary  parameters.  For  unit  twist  per  unit  length  the 
closest  approximation  to  the  stress  function  is  that  for  which  the 
surface  integral 


A 


is  a minimum.  Vhen  <f>  satisfies  this  condition  we  have 

A 


The  Sitz  method  will  now  be  used  in  conjunction  with  a variational 
technique  in  a manner  similar  to  that  used  by  Duncan'.  The  double-wedge 
section  and  the  position  cf  the  origin  and  axes  are  shown  in  Figure  1 . 

In  considering  the  region  O^BB'  it  is  convenient  to  have  the  origin  at 
01t  and  similarly  at  02  for  the  region  OgBB' . A parabolic  variation  of 
the  stress  function  across  the  thickness  of  the  section  is  assumed,  so 
that  in  the  region  C^BB' 

<t>  = (m*  x2  - y^)  & fl  (10) 

and  in  the  region  OgBB1 

<t>  = (*2  4 -y*) s f2  (11) 

In  the  above  equations  f 1 and  f2  are  functions  of  x1  and  Xg  and 
they  will  be  ohosen  to  make  the  surface  integral  K a minimum. 

Substituting  equations  ( 10)  and  (11)  in  equation  (8)  and  integrating 
with  respect  to  y across  the  thickness  gives  K as  the  sum  of  two  integrals 
of  x.j,  f!f  and  Xg,  f2,  f2.  Variations  8f^  in  f^  and  5f2  in  f2  will 
give  rise  to  a variation  6K,  and  for  K to  be  a minimum  6K  must  vanish, 
whence 


/K 
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6K 


■ /’  *i 15  (1-?>  fi  • 10  4 *1 f;  - 2 *5 - 5'  6,i  “i 


* ^ »2  / *2  is  <’-2>  V 10  ■*  *2  V’  »2  4 f"r5<  if2 


+ ~ [{5  m.,  f^)  +t  f;(0l)|  + [5  mg  f2(c2)  +t  f^(o2)|  6f2(c2)] 


(12) 


The  variations  5f and  6f2  are  quite  arbitrary  provided  there  is 
continuity  at  HB' , i.e. 


fl(c.|)  = f2(c2) 

St^Cj)  = 6f2(c2) 


03) 


and  the  expressions  under  the  integral  signs  in  equation  (12)  must 
therefore  vanish.  Similarly  the  expression  in  square  brackets  in  equation 
(12)  must  vanish  subject  to  condition  (13).  The  solution  of  these 
equations  is:- 


'2  = 


Ar,/2m< 


V2m2 
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Substitution  of  equation  (lif)  in  equations  (9),  (10)  and  (H) 
and  integrating  gives 
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apeenux  II 


Calculation  of  upper  limit 

The  method  for  obtaining  an  upper  limit  is  based  on  the  static 
analogue  of  Kelvin’s  theorem:-  "The  strain  energy  of  a structure  corres- 
ponding to  a given  deformation  is  less  than  if  the  freedom  had  been 
limited  by  the  introduction  of  constraints".  The  given  deformation  is 
assumed  to  be  a unit  twist  per  unit  lei^th  and  the  internal  constraints 
are  those  necessary  to  inpose  a chosen  warping  w of  the  cross-section. 

The  position  of  the  oentre  of  twist  is  arbitrary  since  it  may  be 
altered  b£  a rigid  body  movement^,  but  if  it  is  chosen  to  be  at  the 
point  (x,y ) the  strain  energy  per  unit  length  of  cylinder^  is  proportional 
to 


»-//[( <,6> 

A 

and  the  closest  approximation  to  the  warping  function  is  that  for  which 
H is  a minimum.  When  H satisfies  this  cowl  it  ion  we  have 


C 


upper 


GH 


(17) 


The  steps  in  the  analysis  are  similar  to  those  used  in  calculating 
the  lower  limit.  It  is  convenient  to  let  the  section  twist  about  the 
centre  C,  but  in  considering  the  region  0-]BB'  it  is  convenient  to  have 
the  origin  at  0-j,  and  similarly  at  02  for  the  region  OgBB' . A linear 
variation  of  the  warping  funotion  across  the  thickness  of  the  section 


is  assumed,  so  that  in  the  region  O^BB' 

w = yi  e1 

(18) 

and  in  the  region  O^B' 

* = y2  «2 

(19) 

In  the  above  equations  g,  and  g2  are  functions  of  and  x2  and 
they  will  be  ohosen  to  make  the  surface  integral  H a minimum. 

With  the  origins  at  01  and  02  for  the  two  parts  of  the  double 
wedge,  equation  ( 16)  becomes 


/H 
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-C1  ?1x1 


H - / J [yj  (s’  - O2  + (g1  + X1  - ./]  dx1  dyl 


o _inixi 


°2  ®2X2 


/ / ^y2  (g2  " 1)2  + (g2  + X2  ‘ °2)2J  **2  ^2 


o "m2x2 


= -y^  ^ [m2  (g*  - I)2  + 3 x1  (g1  + x1  - Cl)2]  dx1 


2m2  72 


j [m2  x-j  (g^  - I)2  + 3 x2  (g2  + x2  - c2)2]  cbc2 


(20) 


on  integrating  with  reapect  to  y1  and  y2. 

Variations  6g1  in  g1  and  6g2  in  g2  will  give  rise  to  a variation 
6H, and  for  H to  be  a minimum  6H  must  vanish,  whence 


if  ID  1 

6H  = J [3  - 3 c1  + 3 x1  g1  - m2  [3  (g^  - 1)  + x1  glfj]  Sg.jdx., 

o 


3 Xg  c2  + 3 x2  g2 


2 

m2 


(s'z-  ^ + x2  82J]  6«2dx2 


[{S*  (°i)  - ■>]  8S-J  (o,)  + [g^  (c2)  - 1]  6g2  (o2)) 


(21) 


The  variations  &s1  and  6g2  are  quite  arbitrary,  apart  from  contin- 
uity at  BB' , so  that  eaoh  of  the  expressions  in  square  brackets  under 
the  integral  signs  in  equation  (21)  vanish.  The  last  expression  in 
square  brackets  in  equation  (21)  will  vanish  provided  there  is  continuity 
at  BB' , i.e,  provided 


= .-g2(c2) 

Sglto^  = -dg^Og) 


(22) 


The  minus  signs  in  equation  (22)  are  beoause  of  the  reversed  directions 
of  y1  and  y2. 
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The  solution  of  these  options  is 


and 


go 

-2  = 1 
c2 


•(3  a 


(24) 


where 


(m1  + m2)  {m,,-  mg*  p2  (1  - m^)] 
(P1  + P2) ( 1 “ mi)(1  “ 


and 


(®1  + m2^m2  ‘ + Pj  (1  - m1m2^ 

(P1  + PgK1  ' m^)(l  - mp 


Substitution  of  equations  (23)  and  ( 24 ) in  equations  (16)  and  ( 17) 
and  integrating  gives 

Got3  ( 1 ^k(lQaj)  .1,(10  m|) 

upper"  12  \m1  + m2/L  (1  . B2)2  (i  _ m2)2 


+ * m1  m2  (P1  B1  + P2  B2^  ’ 8 m1  m2  (7~f  + TT2^)]  (25) 
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